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ABSTRACT 

With a primary goal of conducting precision weak lensing measurements from space, the 
COSMOS* survey has imaged the largest contiguous area observed by the Hubble Space Telescope 
(HST) to date using the Advanced Camera for Surveys (ACS). This is the first paper in a series 
where we describe our strategy for addressing the various technical challenges in the production 
of weak lensing measurements from the COSMOS data. We first construct a source catalog from 
575 ACS/WFC tiles (1.64 degrees 2 ) sub-sampled at a pixel scale of 0.03". Defects and diffraction 
spikes are carefully removed, leaving a total of 1.2 x 10 6 objects to a limiting magnitude of 
F814W — 26.5. This catalog is made publicly available. Multi-wavelength follow-up observations 
of the COSMOS field provide photometric redshifts for 73% of the source galaxies in the lensing 
catalog. We analyze and discuss the COSMOS redshift distribution and show broad agreement 
with other surveys to z ~ 1. Our next step is to measure the shapes of galaxies and to correct 
them for the distortion induced by the time varying ACS Point Spread Function and for Charge 
Transfer Efficiency effects. Simulated images are used to derive the shear susceptibility factors 
that are necessary in order to transform shape measurements into unbiased shear estimators. For 
every galaxy we derive a shape measurement error and utilize this quantity to extract the intrinsic 
shape noise of the galaxy sample. Interestingly, our results indicate that the intrinsic shape noise 
varies little with either size, magnitude or redshift. Representing a number density of 66 galaxies 
per arcminute 2 , the final COSMOS weak lensing catalog contains 3.9 x 10 5 galaxies with accurate 
shape measurements. The properties of the COSMOS weak lensing catalog described throughout 
this paper will provide key input numbers for the preparation and design of next-generation wide 
field space missions. 

Subject headings: cosmology: observations - gravitational lensing - large-scale structure of Universe 
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1. Introduction 

As we look towards distant galaxies, fluctua- 
tions in the intervening mass distribution cause a 
slight, coherent distortion of their intrinsic shapes. 
This effect, known as weak gravitational lens- 
ing, has been used for more than a decade to 
probe the cosmography and the growth of struc- 
ture (for a review, see iBartelmann fe Schneider 
2001). Although technically challenging because 
the weak lensing signal is minuscule and buried 
in a considerable amount of noise, this field has 
shown substantial progress due to the advent of 
high resolution space based imaging, the prolifer- 
ation of wide-field multi-color surveys, and a de- 
termined effort to improve image analysis meth- 
ods and minimize systematic err ors through the 
Shea r TEsting Program (STEP; iHevmans et al 
2006t iMassey et alJExKih . 



ion: 



Historic ally first observed only around clus- 
ter cores ( Tyson et al. 1990T ). weak lensing has 



emerged as a versatile and effective technique 
to probe th e mas s distribution of clusters (e.g., 



Kneib et al.1 120031 ). to measure the clustering 



straints on the matter density parameter f£ 
and the a mplitude of the matte r power spectrum 



e .g.. iHoekstra et all 120061: ISemboloni et al 



20061: ISchrabback et al. 120061: iHetterscheidt et al 
20061 : Ijarvis et all 20061 ). Most applications cur- 
rently only use the first order deformation in- 
duced by the mass distribution, but novel tech- 
niques are under development to take into ac- 
count second order deformations (also called flex- 



Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by AURA Inc, under NASA 
contract NAS 5-26555; also based on data collected at: the 
Subaru Telescope, which is operated by the National As- 
tronomical Observatory of Japan; the European Southern 
Observatory, Chile; Kitt Peak National Observatory, Cerro 
Tololo Inter-American Observatory, and the National Op- 
tical Astronomy Observatory, which are operated by the 
Association of Universities for Research in Astronomy, Inc. 
(AURA) under cooperative agreement with the National 
Science Foundation; the National Radio Astronomy Obser- 
vatory which is a facility of the National Science Founda- 
tion operated under cooperative agreement by Associated 
Universities, Inc ; and the Canada-France-Hawaii Telescope 
operated by the National Research Council of Canada, the 
Centre National de la Recherche Scientifique de France and 
the University of Hawaii. 



Goldberg fc Baconl 120051: iBacon et all 120061: 



Nataraian et al. 1998: 


Hoekstra et al.ll2004; Sheldon et al. 


2004; Mandelbaum et al. 


2005), and to put con- 



Okura et al.ll2006l ; lGoldberg fc Leonardll2006l ) and 
may prove to be more efficient probes of compact 
structures such as galaxies and groups of galaxies. 
Weak lensing measurements are particularly pow- 
erful when combined with the knowledge of the 
three dimensional galaxy distribution. Sophisti- 
cated lensing 'tomography' techniques that utilize 
redshifts to analyze the 3D shear field are a sen- 
sitive probe of the growth of str ucture and the 
equat ion of state of dark energy dJain fc Taylor 
120031 iBernstein fc Jainl 12004 IBacon et al.ll2005l ). 
Applied in many different ways, weak lensing tech- 
niques unravel the mass distribution of structures 
and their evolution in the universe. 

High quality measurements of weak shear de- 
pend on the accurate determination of the shapes 
and redshifts of distant, faint galaxies. The COS- 
MOS program has imaged the largest contiguous 
area (1.64 degrees 2 ) with the Hubble Space Tele- 
scope (HST) to date using the Advanced Camera 
for Surveys (ACS) Wide Field Channel (WFC). 
The Full Width Half Maximum (FWHM) of the 
Point Spread Function (PSF) of the ACS/WFC 
is 0.12" at the detector, yielding a much bet- 
ter resolution of small galaxies than ground-based 
surveys, which are typically limited by the seeing 
to a PSF of FWHM ~ 1". Shape measurements 
also benefit from ACS/WFC imaging compared to 
ground-based observations because smaller correc- 
tions are required for the PSF and the shear mea- 
surements are less diluted by PSF smearing. The 
imaging quality and unprecedented area of the 
COSMOS ACS/WFC data combined with exten- 
sive follow-up observations at other w avelengths 
(|Scoville et al.ll2007t Ikoekemoerll2007[) to provide 



accur ate photometric redshifts (jMobasher et al 



20071 ). make COSMOS a unique data set for weak 
lensing studies. 

To exploit the weak lensing potential of the 
COSMOS ACS/WFC data, a carefully designed 
catalog of resolved galaxies with shape measure- 
ments must be extracted from the imaging data. 
The challenges and requirements of such a cata- 
log are the following. First, the large survey size 
makes a robust automation of catalog generation 
essential. Second, the lensing sensitivity increases 



1 Before convolution with the detector pixels, the intrinsic 
width of the F814W PSF is 0.085" 
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with the number density of resolved faint galax- 
ies. Thus, it is important to detect all galaxies 
to faint magnitudes while taking care to minimize 
spurious detections that will add noise to weak 
lensing measurements. Third, the high spatial res- 
olution of the ACS/WFC allows for an excellent 
separation of close pairs. Given accurate deblend- 
ing and a high number density of galaxies, one 
can expect to measure shear statistically on sub- 
arcminute scales, where baryonic physics may be- 
gin to influence the dark matter distribution. As 
the COSMOS- ACS/WFC data set is likely to be 
the only large space-based lensing survey until the 
launch of next-generatio n wide field space mis sions 
such as snap dor dune (jRefregier et al.l l2006) . the 
knowledge acquired through the COSMOS data 
will be unique and of crucial importance for the 
preparation and design of these future missions. 

This is the first paper of a series describing the 
galaxy selection, the galaxy shape measurement, 
the lensing analysis, and the cosmological interpre- 
tation of COSMOS data. Details regarding PSF 
corrections as well as tests for systematic effects 
are presented in th e second paper of this series 
( Rhodes et al.l 2007 ) . A three dimensional cosmic 
shear anal ysis is presented in the third paper of 
this series (jMassev et al.l 120071) . Finally, a fourth 
paper presents high resoluti on dark matter mas s 
maps of the COSMOS field dMassev et al.ll2007ah . 

In this paper, we describe our methods for con- 
structing a galaxy catalog from the ACS/WFC 
data, to be used in subsequent weak lensing work 
with COSMOS. Our goal is to produce a catalog 
of galaxies with photometric redshifts and PSF- 
corrected shape measurements, free of contami- 
nating stars, cosmic rays, diffraction spikes, and 
other artifacts. The paper is organized as follows. 
In £j2] we present the data. In $3] we describe the 
pipeline that locates and measures the properties 
of all detected objects. In we assess the qual- 
ity of the data and analyze the COSMOS redshift 
distribution. In Sj5] and $6] we present the PSF 
and CTE correction schemes, the shape and shear 
measurement methods, and our selection criteria 
for the final lensing catalog. In $7] we extract the 
intrinsic shape noise of the galaxy sample as a 
function of redshift and discuss the implications 
for future weak lensing surveys. Where necessary, 



we assume a standard cosmological model with 
il M = 0.3, fl A = 0.7, H = lOO/i km s" 1 Mpc" 1 
and h = 0.7. 

2. The COSMOS ACS Data 



The COSMOS HST ACS field (iScoville et al 



2007; iKoekemoerl 120071 1 is a contiguous 1.64 



2 http:/ /snap, lbl.gov 



degrees 2 , centered at 10:00:28.6, +02:12:21.0 
(J2000). Between October 2003 and June 2005 
(HST cycles 12 and 13), the region was com- 
pletely tiled by 575 adjacent and slightly overlap- 
ping pointings of the ACS/WFC (see Figure [T]). 
Images were taken through the wide F814W filter 
("Broad I"). The camera has a 203" x 203" field of 
view, covered by two 4096x204 8 CCD chips with a 
native pixel scale of 0.05" dFord et al.ll2003l ). The 
median exposure depth across the field is 2028 
seconds (one HST orbit). At each pointing, four 
507 second exposures were taken, each dithered 
by 0.25" in the x direction and 3.08" in the y di- 
rection from the previous position. This strategy 
ensures that the 3" gap between the two chips is 
covered by at least three exposures and facilitates 
the removal of cosmic rays. Pointings were taken 
with two approximately 180° opposed orientation 
angles (PA_V3= 100 ± 10° and 290 ± 10°). In this 
paper we use the "unrotated" images (as opposed 
to North up) to avoid rotating the original frame 
of the PSF. By keeping the images in the default 
unrotated detector frame, they can be stacked to 
map out the observed PSF patterns. For sim- 
ilar reasons, we perform detection in individual 
ACS/WFC tiles instead of on a larger mosaic 
(where the orientation of the PSF frame would 
be unknown). Figure [U shows a COSMOS/ACS 
pointing with the bright detections (F81AW < 23) 
and the masking of stars, asteroid trails and image 
defects (see ^3.5|) . 

To build our catalog, we use version 1.3 of the 
"unrotated" ACS/WFC data which h as been spe- 
cially reduced for lensing purposes fsee lKoekemoer 
2007, for technical details). Image registration, 
geometric distortion, sky subtraction, cosmic ray 
rejection and the final combination of the dithered 
images were performed by the M ultiDrizzle al- 
gorithm (Koekemoer et al.ll2002l ). As described in 
([Rhodes et al.l 120071) . the MultiDrizzle param- 
eters have been chosen for precise galaxy shape 
measurement in the co-added images. In par- 
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Fig. 1 . — Date of observation for each of the sur- 
vey pointings. The PSF of the ACS/WFC varies 
on timescales that are much shorter than the pe- 
riod over which COSMOS was observed. 



ticular, a finer pixel scale of 0.03"/pix was used 
for the final co-added images (7000x7000 pixels), 
even though this implies more strongly correlated 
pixel noise (see H3.8[) . Hereafter, when we refer to 
pixels, we will assume a pixel scale of 0.03"/pix. 
Pixelization acts as a convolution followed by a 
re-sampling and, although current shear measure- 
ment methods can successfully correct for convolu- 
tion, the formalism to properly treat re-sampling 
is still under development for the next genera- 
tion of me thods . Again f ollowi ng the recommen- 
dations of Rhodes et al.l (2007J), a Gaussian and 
isotropic multidrizzle convolution kernel was used, 
with scale=0.6 and pixf rac=0.8, small enough 
to avoid smearing the object unnecessarily while 
large enough to guarantee that the convolution 
dominates the re-sampling. This process is then 
properly corrected by existing shear measurement 
methods. 

The ACS/WFC CCDs also suffer from imper- 
fect charge transfer efficiency (CTE) during read- 
out. As charges are transferred during the read- 
out process, a certain fraction are retained by 



charge traps (created by cosmic ray hits) in the 
pixels. This causes flux to be trailed behind ob- 
jects as the traps gradually release their charge, 
spuriously elongating them in a coherent direc- 
tion that mimics a lensing signal. Since this ef- 
fect is produced by a fixed number of charge traps 
within the CCD substrate, it affects faint sources 
(with a larger fraction of their flux being trailed) 
more than bright ones. This is an insidious effect 
that mimics an increasing shear signal as a func- 
tion of redshift, and prevents the traditional way 
of dealing with the calibration of faint galaxies in a 
lensing analysis by looking at bright stars. Ideally 
this effect would be corrected for on a pixel-by- 
pixel basis in the raw images unfortunately our 
current physical understanding of this effect is in- 
sufficient and a more indepth analysis in still un- 
derway. The CTE effect can be quantified suffi- 
ciently well however that, in a first step, we can 
adopt a post-processing correction scheme based 
on an object's position, flux, and date of observa- 
tion. Fu rther details regardin g this model can be 



3. The COSMOS ACS Galaxy Catalog 

In this section, we discuss the construction of 
the COSMOS ACS/WFC source catalog. This 
catalog is carefully cleaned of defects and arti- 
facts and is made publicly available through the 
Infrared Science Archive (IRS A) database 1 . 

3.1. Detection Strategy 

We use Versio n 2.4.3 of the SExtract or pho- 
tometry package ( Bertin fc Arnoutel 19961 ) to ex- 
tract a source catalog of positions and various 
photometric parameters. In the construction of 
this catalog, our main concern is to pick out the 
small faint objects that contain most of the lens- 
ing signal. The detection strategy that we there- 
fore adopt is to configure SExtractor with very 
low thresholds (even if this leads to more false de- 
tections in the catalog) and to control our sam- 
ple selection via subsequent "lensing cuts" (see 
Sect [6]). We hope to thus reduce unknown selec- 
tion biases introduced by the SExtractor detection 
algorithm. When configured with low detection 
thresholds however, SExtractor also inevitably a) 



1 http://irsa.ipac.caltech.edu/Missions/cosmos.html 
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Fig. 2. — COSMOS pointing acsJJ395836+0141_unrot_sci_12.fits with bright detections, masking and edge 
definition indicated. Adjacent images have a sizable overlap (shown here by the smaller magenta box) which 
allows us to discard detections on the boundaries of each tile (defined by the larger blue box) without losing 
any objects in the final concatenated catalog. The automated masking of the diffraction spikes around bright 
stars (F814 < 23) is pictured here by the red polygons. The basic shape of the star masks is predefined 
and then scaled with the magnitude of the star. The green rectangles correspond to the manual masking 
of aster oid trails and vario us other image defects. Toward the right edge of this image, the dwarf galaxy 
Ll-099 (jlmpev et all [1996) is identified and flagged in a special category. Bright galaxies with F814 < 23 
are depicted by blue ellipses and bright stars with F814 < 23 by red circles. 
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overdeblends low surface brightness spirals and 
patchy irregulars, b) deblends the outer features 
of bright galaxies, c) detects spurious objects in 
the scattered light around bright objects and d) 
undcrdeblends close pairs. 

The correct detection of close pairs enables lens- 
ing measurements on very small scales. However, 
overdeblending and spurious detections adds noise 
to these measurements. In particular, false detec- 
tions around bright objects can have quite high 
signal-to-noise (S/N) values and are not trivial to 
remove with lensing cuts ( see $6]). The method 
presented here is a partial solution to b) c) and 
d). The overdeblending of low surface brightness 
and patchy galaxies remains a difficult problem, 
however, especially for high resolution imaging, 
and calls for an improvement of existing detec- 
tion algorithms. With the advent of high resolu- 
tion multi- wavelength surveys, a possible solution 
forward would be to incorporate color and mor- 
pholo gical information in to the detection process 
(e.g., iLupton et al"1l200ll) . 

While this overdeblending problem persists in 
our catalog, it affects less than 1% of the objects. 
This problem is furthermore mitigated by the cen- 
troiding process during the shape measurement 
stage. Indeed, objects for which the centroid algo- 
rithm fails to converge, which will often be the case 
for overdeblended features, are discarded from the 
catalog. To remedy the remaining problems b) c) 
and d), we adopt and improve the meth od (known 
as the 'Hot- Cold' method) employed in Rix et all 
(2004). In this method, we run SExtractor twice, 
once with a configuration optimized for the detec- 
tion of only the brightest objects ("cold" step) and 
then again with a configuration optimized for the 
faint objects ("hot" step). This double extraction 
helps improve the detection of close pairs. The two 
samples are then merged together to form the final 
catalog and masks are created around the bright 
detections minimizing the effects of c) and d). 

For the "hot" and "cold" steps we vary four 
main parameters to optimize the detection: 1) 
detect_threshold, the minimum signal-to- 
noise per pixel above the background level, 2) 
MIN_AREA, the number of contiguous pixels ex- 
ceeding this threshold, 3) back_SIZE, the mesh 
size of the background map, 4) deblend_nthres 
and DEBLEND_MiNCONT, the parameters regulat- 
ing deblending. In both cases, the data are filtered 



prior to detection by a 5 pixel (0.15") Gaussian 
filtering kernel. Our choice of parameters for both 
steps is provided in Table [TJ 

The two-step method also allows one to adjust 
the estimation of the background map according 
to the typical size of objects one expects to de- 
tect, improving detections with SExtractor. The 
background map is constructed by computing an 
estimator for the local background on a grid of 
mesh size back_size. We adjust back_size so as 
to capture the small scale variations of the back- 
ground noise while keeping it large enough not to 
be affected by the presence of objects. 

For each exposure, a weight map is produced 
by MultiDrizzle describing the combined noise 
properties of the rea d-out, the d a rk cu rrent and 
the sky background ( Koekemoer ( 20071 )). These 
maps describe the noise intensity at each pixel and 
are used to account for the spatial-dependent noise 
pattern in the co-added image with the SExtractor 

WEIGHT_IMAGE option set to WEIGHT MAP. 

Each ACS/WFC pointing consists of four 
slightly offset, dithered exposures making cosmic 
ray rejection more difficult and detection more un- 
reliable on the boundaries of each tile where there 
are fewer than four input exposures. Because ad- 
jacent images overlap sufficiently, we can trim the 
edges of the images without actually removing 
data (see Figure [2]). 

3.2. Bright Object Detection 

In the first step we detect only the brightest 
and largest objects in the image, with 140 or more 
contiguous pixels (corresponding to a diameter of 
0.4" for a circular object) rising more than 2.2 
sigma per pixel above the background level. The 
back_SIZE parameter is set to 12", or 30 times the 
diameter of the smallest objects detected. The de- 
tection threshold and the deblending parameters 
deblend_nthres & deblend_MINCONT are cal- 
ibrated heuristically on several images to separate 
close pairs as much as possible without deblend- 
ing patchy, extended spiral galaxies. Because faint 
objects are captured in a second run, we are free 
to choose the value of the detection threshold dur- 
ing this step. We found that this flexibility greatly 
helped to calibrate the parameters that optimize 
the deblending. Indeed, if detect_threshold 
is set to a low value (for example, to detect faint 
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Table 1 

sextractor configuration parameters 



Parameter 


Bright objects 


Faint objects 


DETECT.MINAREA 


140 


18 


DETECT. THRESH 


2.2* 


1.0 


DEBLEND. NTHRESH 


64 


64 


DEBLEND.MINCONT 


0.04 


0.065 


CLEAN.PARAM 


1.0 


1.0 


BACK.SIZE 


400 


100 


BACK.FILTERSIZE 


5 


3 


BACKPHOTO.TYPE 


local 


local 


BACKPHOTO.THICK 


200 


200 



* Because of correlated noise (sec § !3.8t . the effective thresh- 
old levels are DETECT.THRESH ~ 1.25 for bright objects and 
DETECT.THRESH ~ 0.57 for faint objects. 



object), close pairs will be detected as a single ob- 
ject and are difficult to deblend. Figure [3] illus- 
trates the improvements of this two-step method 
compared to a single-step method. During this 
first step, all pixels associated with a detection 
are recorded by SExtractor in an image called a 
"segmentation map". These segmentation maps 
are used at a later stage to merge the bright and 
the faint catalogs (see N3.4p . This first catalog of 
bright objects is referred to as C co id- 

3.3. Faint Object Detection 

In the second step, we configure SExtractor to 
pick up the small, faint objects, taking care to 
choose the detection parameters to be less con- 
servative than any subsequent lensing cuts (see 
fJ5]). min_AREA is set to 18 pixels (corresponding 
to a diameter of 1.2 times the FWHM of the PSF) 
and the detection threshold is one sigma above 
the background level. As objects detected at this 
step are smaller, the background estimation can 
be improved by refining the mesh size of the back- 
ground map and setting BACK_SIZE to 100 pixels, 
or 20 times the diameter of the smallest objects 
detected. This second catalog is referred to here- 
after as Chot- 

3.4. Merging the Two Samples 

The final catalog is obtained by merging the de- 
tections from C co id and Chot, keeping all objects in 
C co id and only the objects from Chot not detected 
in C co id- To determine which objects to discard 
from Chot, we use the segmentation maps created 
during the bright detection step. To begin with, 




Fig. 3. — This figure illustrates the difficulty 
of correctly deblending close pairs while keeping 
patchy spirals with strong star forming regions in- 
tact. Squares indicate detections from the faint 
step and circles indicate detections from the bright 
step. The top two panels show three objects that 
are not detected or incorrectly deblended by the 
faint step but that are picked by two step method. 
The bottom panel and the arrow towards the star 
in the upper left panel show that even with this 
method, a perfect configuration is still difficult to 
reach. 

we enlarge the flagged areas in these segmentation 
maps by approximately 20 pixels (0.6"). We then 
discard all objects from Chot for which the cen- 
tral pixel lies within a flagged area of these maps. 
Thus, we remove duplicate detections and create a 
mask around all bright objects, immediately clean- 
ing the catalog of a certain number of spurious 
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detections. By a visual inspection of the data, 
we estimate that this method solves about half of 
the deblending problems that we observe (exclud- 
ing the low surface brightness galaxies). The final 
catalog of raw SExtractor detections C\ contains 
1.8 million objects in total (see Table [2]). 

3.5. Cleaning the Catalog 

Great care was taken to mask unreliable regions 
within images and to remove false detections from 
the catalog, especially those that can mimic a lens- 
ing signal. As illustrated in Figure [2j an auto- 
matic algorithm was developed to define polygo- 
nal shaped masks around stars with F814W < 19 
(the limit at which stars saturate in the COSMOS 
images), with a size scaled by the magnitude of the 
star. Objects near bright stars or saturated pix- 
els were masked to avoid shape biases due to any 
background gradient. All the images were then vi- 
sually inspected. In a few cases the automatic al- 
gorithm failed (very saturated stars for which the 
centroid of the star is widely offset) and the stel- 
lar masks were corrected by hand. Other contam- 
inated regions of the images were also masked out, 
including reflection ghosts, asteroids, and satel- 
lite trails. Astronomical sources such as HII re- 
gions around bright galaxies, stellar clusters, and 
nearby dwarf galaxies, were also flagged and re- 
moved from the lensing catalog. 

Objects with double entries in the catalog (from 
the overlap between adjacent images) are identi- 
fied and the counterpart with the highest SExtrac- 
tor flag (indicating a poor detection) is discarded, 
leaving a catalog of unique objects. However, the 
duplicated objects from overlapping regions are a 
valuable asset for consistency checks and are used, 
for example, to check the galaxy shape measure- 
ment error (see $7] and Figure [T7|). 

The final clean catalog (C2) is free of spurious 
or duplicate detections and contains 1.2 million 
sources in total (see Table [2j . 

3.6. Star-Galaxy Classification 

The correct identification of stars has two im- 
plications for the lensing analysis. First, bright 
stars are useful for PSF modelling and second, 
stars must be correctly identified in order to ap- 
ply our automatic masking algorithm of diffraction 
spikes. A robust star-galaxy classification is thus 



necessary. 

SExtractor produces a continuous stellar 
classification index parameter ranging from (ex- 
tended sources) to 1 (point sources). This index 
has two drawbacks: first the definition of the di- 
viding line is ambiguous and second, the neural- 
network classifier used by SExtractor is trained 
with ground-based images and is therefore only 
valid for a sample of profiles similar to the orig- 
inal training set. With space-based images, this 
index becomes difficult to interpret, as illustrated 
in Figure [4] which depicts our star selection (de- 
scribed below) within the CLASS_STAR/mag_AUTO 
plane. 

We therefore test two alternative methods to 
classify point sources and galaxies, one based on 
the SExtractor parameter mu_max (peak sur- 
face brightness above the background level) and 
the ot her based on the half-light radius , Rhl 
(e-g-, 



Peterson et all 1 19791: iBardeau et ail (20051) . 



Both methods are motivated by the fact that the 
light distribution of a point source scales with 
magnitude. Point sources therefore occupy a 
well-defined locus in a mu_max/mag_AUTO or 
a Rhl/mag_AUTO plane. Figure [5] shows how 
we can use this property to define stars (ID=2) 
and galaxies (ID=1) reliably up to a magnitude 
of F8141F ~ 25. At fainter levels, the clas- 
sification begins to break down and the point 
sources become indistinguishable from the small 
galaxies. We find that the two methods agree 
very well, within 1% at magnitudes less than 
^8141^ = 24 and within 2% at magnitudes less 
than F8141F = 25. The small difference arises 
mainly from a misclassification of objects by the 
Rhl method because of the presence of a close pair 
that distorts the estimation of Rhl. Overall, the 
MU.MAX method proved to be more robust and has 
the advantage of a tighter correlation of the stellar 
locus and a clear break indicating the magnitude 
at which the stars saturate (mag_AUTO ~ 19). 
Moreover, with this method, a surface brightness 
cut at the faint end of the stellar sample is trivial 
to implement (at faint magnitudes, the catalog is 
surface brightness limited). The performance of 
this star galaxy separation scheme will be ana- 
lyzed in more detail in § 14.11 

Using the MU_MAX method, we also define a set 
of objects that are more sharply peaked than the 
PSF, which is obviously non-physical. A visual 
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Table 2 

sumary of the construction of the acs lensing catalog 



Catalog of raw SExtractor detections: Ci 


Number 


Percent of Ci 


Total number of objects in C\ 
Number of Hot (faint) detections 
Number of Cold (bright) detections 


1.8 X10 6 
1.6x10 s 
2.2x10 s 


88% 
12% 


Details of the cleaning process 


Number 


Percent of C\ 


Number of objects within the noisy border of a tile 

Number of Hot detections with central pixel in Cold segmentation map 

Number of objects within automatically defined star masks 

Number of objects within manually defined masks 

Number of objects detected more than once in adjacent tiles 


3.2 x10 s 
2.0x10 s 
2.4xl0 4 
4.1X10 4 
6.6X10 4 


17% 
11% 
1% 
2% 

4% 


Catalog cleaned of image defects: C2 


Number 


Percent of C2 


Total number of objects in C2 
Number of galaxies (ID — 1) 
Number of point sources (ID — 2) 
Number of fake detections (ID — 3) 


1.2x10 s 
1.1x10 s 
2.8X10 4 
1.7 xlO 4 


96 % 
2 % 
2% 


ACS galaxies from C 2 with FSUWab < 26.5: C 3 


Number 


Percent of C3 


Total number of objects in C3 

Number of galaxies with a counterpart in the photometric catalog 

Number of galaxies that have been matched but that arc in ground based masks 

Number of galaxies for which the redshift code did not converge 

Total number of galaxies with accurate photometric rcdshifts, C3 


7.0x10 s 
6.0x10 s 
8.3X10 4 
l.lxlO 4 
5.0 x10 s 


85% 
12% 
1.7% 
71% 


The final COSMOS ACS lensing catalog: C 4 


Number 


Galaxy number density 


Total number of galaxies in final lensing catalog 
Number of galaxies with accurate photometric rcdshifts 


3.9x10 s 
2.8x10 s 


66 arcmin 2 
48 arcmin 2 
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Fig. 4. — The SExtractor stellar index 
(CLASS_STAR) for our point source selection based 
on the peak surface brightness of objects. Grey 
points show the corresponding galaxy sample. A 
point source selection of the form CLASS_STAR > 
0.8 for example, will miss a certain number of 
bright stars (F814W < 17) and will grossly mis- 
classify compact galaxies at 22 < F81AW < 25. 



Fig. 5. — Classification of point sources, galax- 
ies and artifacts within the mu_max/ mag_AUTO 
plane. Point sources follow the PSF and are de- 
limited by the solid region. Objects that are more 
sharply peaked than the PSF are contained in the 
dashed region and are considered to be artifacts. 
For clarity only a 2% random selection of all ob- 
jects are included in this plot. 



inspection finds that these objects are mainly ar- 
tifacts, hot pixels, and residual cosmic rays. We 
flag these spurious objects in our catalog (ID=3) 
and remove them for the lensing analysis. 

Averaging over the COSMOS field, we find ~ 15 
stars per pointing with 19 < MAG_AUTO < 23. 
This is an insufficient number to model the PSF 
in individual images using standard interpolation 
techniques. However, it is a sufficient number to 
identify the PSF pattern of each exposure given 
a finite set o f recu rring patterns (see § 15.11 and 
Rhodes et~al"1 (|2007t) for further details). 



3.7. Photometric Redshifts 

In addition to the ACS/WFC (F8UW) imag- 
ing, the COSMOS field has been imaged with Sub- 
aru Suprime-Cam (Bj,Vj,g+,r+,i+, z+,NB816), 
the Canada-French Hawaii Telescope (CFHT) 
(u*,i*) and the KPNO/CTIO (Ks). Details of 
the ground-based obser vations and t h e dat a re- 
duction are presented in ICapak et alj ( 2007 ) and 



XMM-Newton and in the radio with VLA, CSO 
and IRAM. Yet more observations are underway 
including intermediate and narrow-band imaging 
with Subaru Suprime-Cam, deep Infrared imag- 
ing covering 1.0-2.2 microns (WIRCam/CFHT, 
WFCAM/UKIRT and ULBcam/UH2.2), and ob- 
servations with space-based facilities including 
Chandra and Spitzer Space Telescopes. This ex- 
tensive multi-wavelength data-set is a key com- 
ponent to COSMOS weak lensing measurements 
because it allows us to accurately measure the 
COSMOS redshift distribution, to separate fore- 
ground and background structures, and to remove 
contamination from i ntrins ic galaxy alignments 
(jHevmans fc Heavens 20031) and shear-ellipticity 
correlations ([King 2005J). 

Photometric redshifts were determined by 
the COSMOS photometric redshift code with a 
Bayesian prior based on luminosi ty functions and 



Taniguchi et al.l (|2007). Other observations were 



taken in the UV with GALEX, in the X-ray with 



allow ing for internal extinction (jMobasher et al 
120071) . For each galaxy, the entire probability 
distribution P(z), the most likely redshift and 
a confidence level for that redshift is calculated. 
The knowledge of the full P(z) allows us to apply 
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weight to galaxies in weak lensing measurements 
according to the uncertainty in their measured 
redshift. The accuracy of the photometric red- 
shifts are estimated based on extensive simula- 
tions and by comparison with a sample of 958 
galaxies with spectroscopic redshifts measured 
by the ESO/VIMOS instrument as p art as the 
zCOSMOS program ((Lilly et al.ll2007l ). Both the 
luminosity prior and the extinction corrections 
have been shown to improve the accuracy of the 
photometric redshifts when compared to the spec- 
troscopic sample. The accuracy of the current 
photometric redshifts down to F814IU ~ 22.5 is: 



<W(l + «») =0.031 



(1) 



with r\ = 1.0% of catastrophic errors, defined as 
Az/(l + z s ) > 0.15. This relati on scales wi t h mag - 



nitude in a similar fashion to IWolf et alJ (|2004l ). 



The accuracy of the photometric redshifts will con- 
tinue to improve as more data become available 
(in particular the deeper u*, J, K, Spitzer and 
Subaru narrow band data). 

The COSMOS o ptical and near infrared catalog 
(|Capak et al J 120071 ) provides multi-band photom- 
etry for 89% of the COSMOS ACS/WFC galaxies. 
As demonstrated by Figure [BJ the remaining 11% 
of the galaxies for which we lack multi-wavelength 
information (and therefore photometric redshifts), 
are the small galaxies that cannot be detected by 
the ground based imaging. In effect, although the 
SUBARU data is deeper for sources 1" in diam- 
eter, the ACS/WFC imaging will do a better job 
at detecting anything smaller. Because we apply 
a size cut to the the final lensing catalog how- 
ever (see § H]) , many of these small galaxies will be 
discarded from the final analysis. In total, after 
removing the galaxies with unreliable multi-band 
photometry (because they are masked out in the 
ground-based data) as well as those for which the 
photometric redshift code failed to converge, 76% 
of the galaxies in the COSMOS ACS/WFC lens- 
ing catalog (F814W < 26) have photometric red- 
shifts (See Figure [6J. For further discussions on 
the photometric redshifts and the COSMOS red- 
shift distribution, see § 14.31 



3.8. Noise Properties 

The drizzling process introduces pattern-dependent 
correlations between neighboring pixels and can 
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Fig. 6. — In the upper panel, the number counts 
of galaxies that have been correctly matched to 
the ground based catalog (dotted line) are com- 
pared to the total number counts (solid line) . The 
dashed line indicates galaxies for which we con- 
sider the photometric redshift to be reliable. The 
difference between these two curves is primarily 
due to larger masked areas in the ground-based 
data than in the ACS imaging. In the lower 
panel, the sizes of galaxies within the ACS catalog 
(solid line) are compared to the sizes of those that 
have been matched with the ground based catalog 
(dashed line). The objects for which we do not 
have multi-band photometry are small galaxies 
that are detected with ACS but not with ground 
based imaging (seeing ~ 1"). The vertical dashed 
line shows the approximate size cut that we make 
in the final lensing catalog. 

artificially reduce the noise levels in co-added im- 
ages. Noise and error estimates derived from driz- 
zled images will thus tend to underestimate the 
true noise levels of the image. One should in the- 
ory take into account the exact covariance matrix 
of the noise in order to derive error estimates for 
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drizzled images. For our purposes, however, a 
simple scaling of the noise level is each pixel by 
the same constant factor is sufficient. The scal- 
ing factor that we adopt, Fa, has been derived 
for Mu ltiDrizzled images by ICasertano et al. 
(l2000h . In principle Fa is size dependent but 



converges rapidly with increasing size toward an 
asymptotic value given by: 



(«/p)(l - */(3p)) (s<P) 
(l-p/(3«)) (p<s) 



(2) 



where p and s are respectively the pixf rac and the 
scale configuration parameters of MultiDrizzle. 
For the purpose of this paper, we assume Fa to 
be constant and equal to Fa ~ 0.316 (p = 0.8 and 
s = 0.6). By assuming a constant corrective factor 
regardless of size, we make less than a 10% error 
on the noise estimation of the smallest objects de- 
tected. 

As implemented by SExtractor, the formulas 
for the flux and magnitude uncertainties (for both 
auto and the ISO quantities) are given by: 



flux_err = y/ Act 2 + F/g 



MAG_ERR = 



2.5 FLUX_ERR 



In 10 



F 



(3) 



(4) 



where A is the area (in pixels) over which the flux 
F (in ADU) is summed and g is the detector gain. 
To correct the magnitude and flux errors reported 
by SExtractor for the correlated noise, we replace 
<t, the standard deviation of the noise (in ADU) 
estimated by SExtractor, by a I\JFa within the 
above equations. The significance of a COSMOS 
detection, after this correction is applied, is de- 
fined as s/n = flux_auto/fluxerr_auto (see 
Figure [21). 

4. Quality Assessment of the ACS Catalog 
4.1. Galaxy and Stellar Counts 

The same mu_max parameter used to classify 
stars and galaxies can also be used as an indica- 
tion of the background level and the depth of the 
data. For each image i we calculate the mode nii 
of the MU.MAX parameter. We then divide the m s 
into two bins according to the angle of the tele- 
scope with the Sun at the moment of the pointing. 




Fig. 7.— The significance of COSMOS de- 
tections defined as flux_auto/fluxerr_auto 
where fluxerr_AUTO has been corrected for cor- 
related noise. 

The histogram of m, (Figure [8]) for these two bins 
reveals that the depth of the data is bimodal de- 
pending on whether the angle of the sun is less or 
greater than a critical angle of 70 deg. This is vi- 
sually evident when we inspect the density map of 
the very faint objects (Figure [8]). 96 pointings out 
of 575 have a Sun angle less than the critical value 
making them slightly shallower than the average. 

The number counts serve as a check of the ap- 
proximate photometric calibration and the depth 
of the data. Figure [9] shows the counts for galaxies 
and sta rs compared to the r eported HDF F814W 
counts ( Williams et al.lll996h . The magnitudes are 
given in the AB-system. Stars have been sub- 
tracted from the galaxy counts up to F8UW = 
25. At fainter magnitudes their contribution is 
negligible. We plot raw number counts only, i.e. 
we do not correct for incompleteness at the faint 
end. To facilitate comparisons with other surveys, 
we fit the galaxy counts between ^81414^ = 20 
and F814W — 26, to an exponential of the form 
N = B x io Axma 9 where N has units of number 
degree -2 0.5 mag -1 . For the deeper set of images 
(sun angle > 70 degrees), we find A = 0.332 and 
logio(B) = —3.543. The raw galaxy and stellar 
number counts are provided in Table [3) 

We also fit the stellar counts to models as shown 
in Figure [9] The star count predictions have been 
done using the Besangon model of the Galaxy 
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Fig. 8. — The top panel shows the histogram of mi 
(mode of the MU_MAX parameter for each image) 
demonstrating that images with a sun angle less 
than 70 degrees are not quite as deep as images 
with a sun angle greater than 70 degrees. The 
bottom panel shows the density of faint objects 
(27 < F8UW < 28) within the COSMOS field. 
The red circles indicate the pointings for which 
the angle with the sun is less than 70 degrees. 



dRobin et al.ll2003l 2004) and are described in de- 



tail in (jRobin et al.l 120071 ). By extrapolating the 



stellar number counts, we estimate that the galaxy 
catalog has less than a 3% contamination from 
stars at magnitudes greater than 25. The fit to 
stellar models is excellent between F814W = 20 



and FS14W = 25, and at magnitudes less than 
19, we visually inspect the catalog to check that 
the star selection is correct to within 0.5%. This 
small error arises mainly from false detections by 
SExtractor of the diffraction spikes of bright, sat- 
urated stars. 
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Fig. 9. — The top panel shows the galaxy and stel- 
lar number counts as compared to the HDF. The 
dashed curve corresponds to images with a sun- 
angle of less than 70 degrees and the solid curve 
corresponds to images with a sun-angle greater 
than 70 degrees. Poisson error bars are also in- 
dicated but are very small. The bottom panel 
shows the point source selection for the catalog 
compared to stellar models computed from evolu- 
tionary tracks and constrained by local Hipparcos 
data. 



4.2. Completeness 

The probability that a galaxy enters our cat- 
alog will depend on its size and surface bright- 
ness profile. To quantify the completeness and 
detection limits of our SExtractor configuration, 
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Table 3 

Cosmos F814W galaxy and stellar number counts 



FS14W 


Galaxy density 
l°9lo( n ) deg -2 0.5 mag -1 


Stellar density 
l°9lo( n ) deg -2 0.5 mag -1 


20.25 


3.138 ± 1.503 


2.807 ± 1.297 


20.75 


3.323 ± 1.600 


2.865 ± 1.329 


21.25 


3.514 ± 1.691 


2.930 ± 1.361 


21.75 


3.686 ± 1.777 


2.936 ± 1.365 


22.25 


3.853 ± 1.860 


2.981 ± 1.385 


22.75 


4.022 ± 1.945 


3.003 ± 1.401 


23.25 


4.180 ± 2.023 


3.043 ± 1.414 


23.75 


4.352 ± 2.110 


3.081 ± 1.437 


24.25 


4.523 ± 2.196 


3.143 ± 1.463 


24.75 


4.682 ± 2.275 


3.206 ± 1.493 


25.25 


4.814 ± 2.340 




25.75 


4.956 ± 2.412 





Note. — Galaxy counts are derived for the 479 images with a sunanglc 
greater than 70 degrees. Magnitudes arc the SExtractor MAG-AUTO 



we insert fake objects with a Gaussian profile of 
varying FWHM and total magnitude into empty 
regions of an ACS image and test how well these 
objects can be recovered with our pipeline. Each 
artificial source was considered to be correctly de- 
tected if its centroid was within 10 pixels and the 
MAG_AUTO parameter was within 0.5 mag of the 
input value. From this analysis, we determine 
our completeness as a function of magnitude and 
FWHM and the results are shown in Figure [TO] 
The completeness is about 90% for objects with 
a FWHM of 0.2" at F814W=26.6. These values 
should only be used as a rough estimate however, 
as we do not actually model galaxies, but use a 
simple Gaussian profile for artificial objects. 

4.3. The COSMOS Redshift Distribution 

The estimation of the redshifts of galaxies is 
the major astrophysical uncertainty inherent to 
weak lensing methods. To first approximation, 
cosmic shear and tomography are mainly sensi- 
tive to the median redshift of the sources while 
galaxy-galaxy lensing benefits greatly from the 
knowledge of precise spectroscopic redshifts fo r 
the foreground lenses (jKleinheinrich et al. l2005h . 
With the depth and area coverage of COSMOS, 
we surpass the current capability for complete 
spectroscopic follow-up. For most forthcoming 
weak lensing surveys, this will also be the case, 
hence the importance of the photometric red- 
shift technique to measure redshifts for a major- 



ity of the galaxie s, to and beyond today's spec- 



troscopic limits (jllbert et al. 2006). 



COSMOS 

presents a unique advantage, in terms of cur- 
rent weak lensing surveys, of a prodigious multi- 
wavelength follow-up combined with the planned 
measurement of ~ 50000 spectroscop ic redshifts 
by th e ongoing zCOSMOS program (Lilly et al 
20071 ). Upon completion, this data set will pro- 



vide the COSMOS lensing catalog with accurate 
photometric redshifts and will be vital for refin- 
ing and improving the photometric technique in 
preparation for forthcoming weak lensing surveys. 
We present here a first analysis of the COSMOS 
redshift distribution. A more detailed study of the 
systematic trends in the photometric redshifts and 
of their effects on the redshift distribution is be- 
yond the scope of this paper and will be addressed 
elsewhere when more data becomes available. For 
the purposes of this paper, we adopt the mag- 
nitude dependent parametrization of the redshift 
distribution, comm on to many other weak len sing 
studies as given bv lBaueh fc Efstathioul (fl993h. 



n(z,mag) oc z 2 exp[ 
z {mag) = \ 



z (mag) 



1.5] 



(5) 



where z m is the median redshift of the survey as 
a function of magnitude. We calculate z m for the 
COSMOS ACS data by bins of A.F8141F = 0.25 
for 20 < F8UW < 24 and derive the best linear 
fit to z m , given by: 
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Fig. 10.— Completeness of the COSMOS F814W 
catalog as a function of total magnitude and 
FWHM determined by inserting fake objects into 
an ACS image. The thick contours show the per- 
centage of fake objects recovered by SExtractor. 
The thin contours are the lines of constant sur- 
face brightness, in units of mag arcsec -2 , assum- 
ing a Gaussian profile. The grey points represent 
a random sample of objects from the COSMOS 
catalog plotted as a function of MAG_AUTO and 
FWHMJMAGE. The dashed horizontal line indi- 
cates the size of the ACS PSF. Note that the sim- 
ulations only consider objects with Gaussian pro- 
files whereas in reality, the COSMOS objects ex- 
hibit a wide variety of profiles. 



z m = (0.18 ± 0.01) x F8UW - (3.3 ± 0.2) (6) 

The majority of the galaxies for which we have 
no redshift estimate at F814W < 24 are those 
in masked regions and their exclusion from this 
derivation do not affect this results. At F8UW < 
24, our redshift incompleteness is less than 4% and 
the dominant source of error is the photometric 
redshift uncertainty expressed previously in Equa- 
tion!]] We note that a significant number of galax- 
ies fainter then this limit do have photometric red- 
shifts. The fact that these galaxies represent a sta- 
tistically incomplete sample, does not matter for 
some applications. For example, these additional 
ga laxies are used in the cosmic shear measurement 
bv lMassev et all (|2007t l. 




23 
F814W AB 



Fig. 11.— Median redshift of COSMOS compared 
to various photometric redshift surveys. For clar- 
ity, error bars are only shown for the COSMOS, 
CFHTLS and UDF surveys 



In Figure 111! we compare the median redshift 



of COSMOS to the U DF survey dCoe et al.l l20"(ih. 
the CFHTLS su rvey (lllbert et all l20odT the H- 
HDF-N survey (|Capak et all 12004 and SSA22 
(Hu et al. 2004, Capak et al. 2004). Table H 
is a summary of the data and the methods em- 
ployed by these different photometric redshift sur- 
veys. All photometric redshifts have been com- 
puted with a Bayesian prior based on luminosity 
functions. The agreement that we see between the 
various surveys at z < 1 is quite remarkable. At 
z > 1 however, the scatter in Figure fTTI indicates 
the limits of current photometric techniques. Fur- 
ther simulations are clearly necessary in order to 
understand the biases introduced in the redshift 
distribution at z > 1. Although we reserve a full 
discussion for a future paper, we can already high- 
light some of the issues at hand. 

The photometric redshift technique relies on de- 
tecting and measuring the strength of broad spec- 
tral features. These same features are used by 
color selection techniques to select objects at spe- 
cific redshifts. The key features are the 4000A 
break, the Lyman break at 912A;- Lyman absorp- 
tion at I2I6A7 and coronal line absorption be- 
tween 1500-2500A (see Adelburger et. al. 2005). 
Photometric redshifts are very robust if one or 
more of these features are detectable in the avail- 
able data. However, at faint magnitudes the pho- 
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Table 4 

Present date photometric redshift surveys 



Survey 


Area 


Imaging data 


Calibration spectra 


Technique 


COSMOS 


1.67 deg 2 


B i i v i i 9+, r+, i+, z+, JVB816, u*,i*, Ks 


958 


COSMOS/BPZ 


UDF 


11.97 arcmin 2 


B,V,i',z',J* H h 


76 


BPZ 


CFHTLS 


3.2 deg 2 


«*, p', r', i\ z\ J, K 


2867 


Lc Pharc 


H-HDF-N 


0.2 deg 2 


1 .11 A .li. l ./..UK 


2149 


BPZ 


SSA22 


0.2 deg 2 


U- ,B j ,V j ,R c ,I c ,z+, J, H, K, HK' 


452 


BPZ 



a 5.76 arcscc 2 
b 160 arcmin 2 



tometric errors and detection limits are often too 
large to constrain these features. 

For example, problems arise for COSMOS be- 
tween 1.5 < z < 3.2 where the measurable fea- 
tures, the 4000A break, and the coronal line ab- 
sorption features, are difficult to detect. Indeed, 
at these redshifts, the 4000A break is well into the 
IR where it is difficult to obtain deep data. A 
typical object at z ~ 2 will be ~ 1.4 magnitudes 
fainter at I than K. This means objects fainter 
than .F814W > 23.5 are not constrained by the 
present if-band data. At the other end of the 
spectrum the coronal absorption feature has a typ- 
ical strength of ~ 0.15 magnitudes, which requires 
a 25ct detection in u* to accurately differentiate 
from a similar break in z < 0.5 galaxies. With 
the present u* data, this corresponds to objects 
brighter than F814W < 24.5 for typical galaxies. 

In conclusion, the high redshift tail of the COS- 
MOS redshift distribution will be more accurately 
determined with the forthcoming NIR data and 
the future deep zCOSMOS spectroscopy which 
specifically targets the 1.5 < z < 3 region and will 
allow proper calibration down to F814W ~ 24. A 
more detailed analysis of the photometric redshifts 
will be conducted once the new data is available. 

5. PSF Correction and Shear Measure- 
ment 

In this section, we measure the shapes of galax- 
ies and correct them for the convolution with the 
telescope's PSF and for other instrumental effects. 
For each galaxy, we construct an unbiased local 
estimator of the shear and derive the associated 
measurement error. 



5.1. PSF Modelling 

The ACS/WFC PSF is not as stable as one 
might naive ly hope from a spa ce-based camera. 
As shown in lRhodes et al. ( 2007 ). gradual changes 
to both the size and the ellipticity pattern of the 
PSF due to telescope "breathing" , causes the PSF 
to change considerably on timescales of weeks. 
The long period of time over which the COSMOS 
field was observed forces us to take account of 
these variations (see Figure [T]) . Although other 
strategies have been demonstrated successfully for 
observations conducted on a shorter time span, 
it would be ina p propr iate for us to assum e, like 



Lombardi etaD (|2005h or Ijee et all d2005h. that 



the P SF is constant or even, like iHevmans et al 
(2005), that the focus is piecewise constant. 



Fortunately, most of the PSF variation can be 
ascribed to a single physical parameter. Thermal 
expansions and contractions of HST alter the dis- 
tance between the primary and secondary mirrors. 
As the "effective focus" deviates from nominal, the 
PSF becomes larger and more elliptical, with the 
direction of elongation depending upon th e po- 
sition above or below nominal focus (c./. iKristl 
2005). The thermal load on HST is constantly 
changing, in a complicated way, as it passes in 
and out of the shadow of the Earth and is rotated 
to different pointings. 



As described in iRhodes et al.l (|2007l ). we have 
modified version 6.3 of the TinyTim ray-tracing 
program (Krist and Hook 2004) to create a grid 
of model PSF images, at varying focus offsets. By 
comparing the ellipticity of ~ 20 stars in each im- 
age to these, we can determine the image's effec- 
tive focus. Tests of this algorithm on ACS/WFC 
images of dense stellar fields confirm that the best- 
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5.2. Galaxy Shape Measurement 



150.6 150.4 



150.2 150.0 
RA 
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Fig. 12. — Adopted PSF model across the sur- 
vey. The colors correspond to deviations in the 
apparent focus of the telescope away from nom- 
inal (fj,m). These are caused by expansion and 
contraction of the HST due to thermal variations. 
Note that the focus values are clustered and not 
randomly distributed. See Rhodes et al. 2007 for 
more details about the PSF pattern at particular 
focus positions. 



fit effective focus can be repeatably determined 
from a random sample of ten stars brighter than 
F814W — 23 with an rms error less than l/im. 
The effective focus of the COSMOS images are 
shown in Figure [TH An alternative correction 
scheme based on PSF models constructed from 
dense stellar fields has also been suggested by 
~ (2Q03). 



Schrabback et al 



Once images have been grouped by their effec- 
tive focus position, we can combine the few stars 
in each image into one large catalog. We interpo- 
late the PSF model parameters using a polynomial 
fit (of order 3 x 2 x 2 in each CCD s eparately), in 
the u sual weak tensing fa s hion (c.f. iMassey et al. 
20021 ). See iRhodes et al.l (j2007l ) for more details 



We use the shape measuremen t method de- 
velope d for space-based imaging bv lRhodes et al 
(|2000L hereafter RRG). The RRG method has 
been optimized for space-based images with 
small PSFs and has previously been used on 
weak lensing analyses of WFPC2 and STIS 
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data (Rhodes et al. 2001 : iRefregier et al. 
Rhodes et al.l 12004 ). In a manner similar to 
the common "KSB" method (jKaiser et al.lll995l ). 
RRG measures the second and fourth order 
Gaussian-weighted moments of each galaxy: 



Hjkl 



^2 wlXiXj 

Y J wIx l x j x k xi 



(7) 



(8) 



The sum is over all pixels, w is the size of the 
Gaussian weight function, I is the pixel inten- 
sity, and the Xi coordinates are measured in pixels. 
The Gaussian weight function is necessary to sup- 
press divergent sky noise contributions in the mea- 
surement of the quadripole moments. The RRG 
method is well-suited to the small, diffraction- 
limited PSF obtained from space, because it de- 
creases the noise on the shear estimators by cor- 
recting each moment for the PSF linearly, and only 
dividing them to form an ellipticity at the last pos- 
sible moment. 

After the moments have been corrected for the 
PSF, an ellipticity e = (ex, €2) and size measure, 
d, are calculated for each galaxy: 



('1 



('2 



L !J)J 
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{I xx + lyy) 



(9) 
(10) 

(11) 



concerning the PSF modelling scheme. 



Note that the d parameter is a measure of 
galaxy size but that its value will depend on the 
choice of the width of the Gaussian weight func- 
tion, w. 
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5.3. Shear Measurement 



galaxy as a function of S/N. 



The estimator e — (ei,e2) is not yet a shear 
estimator, because it does not respond linearly to 
changes in shear. It must first be normalized by 
a shear susceptibility factor (also known as the 
"shear polarizability" ) , 



7 = {7i>72} 



£ 

G 



(12) 



where the shear susceptibility factor, G, is mea- 
sured from moments of the global distribution of 
£ and other, higher order shape parameters (see 
equation 28 in lRhodes et~ai1l2000h . The RRG for- 
malism does not allow for G to be calculated for 
any individual galaxy. However, G can be cal- 
culated for an ensemble of galaxies by averaging 
over a population's shape moments. Previous in- 
carnations of RRG that have been used to measure 
cosmic shear (Rhodes et al.ll2001t iRefregier et al 



2002; IRhodes et al.ll2004r ) , have made use of a sin- 
gle value of G for the entire survey. Adopting 
this approach for the COSMOS data would yield 
a value of G = 1.13. However, the Shear TEst- 



ing Program (M assev et al.l 120071 ) showed that G 
can vary significantly as a function of object flux, 
whether this be due to evolution in galaxy mor- 
phologies as a function of redshift, or noise in the 
wings of faint galaxies that simply impedes the 
measurement of their radial profiles and higher or- 
der moments. An increase in shear susceptibility 
with obje ct S/N has also bee n seen in KSB-type 
analyses ([Massev et al.ll2004al ). 

Our tests have confirmed that a constant value 
would be insufficiently precise for a survey the 
size of COSMOS, and would particularly affect the 
kind of 3D analysis for which COSMOS is so well- 
suited. We therefore calculate G from the COS- 
MOS data in bins of S/N (see Figure [13]). Sim- 
ulated images of the COSMOS dat a are created 



using th e shapelets-based method of lMassev et al 
(|2004bl ) (see gO]). The variations of G as a func- 



tion of S/N are apparent in the COSMOS data 
are well reproduced by the simulated data. For 
COSMOS galaxies, we find that variations in G as 
a function of S/N are well-fit by 



a = .1.12:, tui4i.ivsi.il ( S//N - 17 ) . (13) 



Adopting the above model, we derive G for each 




Fig. 13. — Interpolation of the shear susceptibility 
factor G. The solid circles show G, calculated in 
bins of S/N, for the COSMOS data. The open 
circles show the same for the simulated COSMOS 
images. The solid line shows the shear suscepti- 
bility model adopted for the data. 

5.4. Calibration via Simulated Images 



Usin g the shapelets-based method of lMassev et al 
(|2004bT ). we have created simulated images with 
the same depth, noise properties, PSF, and galaxy 
morphology distribution as the real COSMOS 
data. A known shear signal was applied to the 
images, which we have then attempted to mea- 
sure using the same pipeline as the data. This 
exercise is similar to the Shear TEsting Program 
(STEP; iHevmans et~ai1l200d ; iMassev et al"1l2006l ) 
but tailored exclusively to COSMOS. 

The simulated COSMOS images are each 4' x 4', 
and contain ~ 500 galaxies after applying the 
same catalog cuts that were applied to the real 
data (see Simulated galaxy morphologies 

are bas ed on those observed in th e Hubble Deep 
Fields ( Williams et al.lfl99^ . Il998f ). parametrized 
as shapelets and randomly rotated/flipped before 
being sheared. Different input galaxies were used 
in each simulated image, as if they were pointing to 
different patches of the sky, to keep them indepen- 
dent. To simplify later analysis, all of the galaxies 
within an image were sheared by the same amount. 
A total of 41 images were made, with shears ap- 
plied in integer steps from —10% to +10% in the 
71 component (while 72 was fixed at zero) and 
similarly for the 72 component. The images were 
then convolved with a model ACS PSF. Again to 
simplify the analysis, this was a constant PSF ob- 
tained from TinyTim. Its (ei,e 2 ) ellipticity is 
(-0.21%, -2.07%) ± (0.14,0.10). No stars were 
included in the simulated images; a separate star 
field was created, from which the PSF moments 
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could be measured. Noise was added to all of these 
images, to the same depth as the COSMOS ob- 
servations, and with a similar (but isotropic) cor- 
relation between adjacent pixels to mimic the ef- 
fects of MultiDrizzle and unresolved background 
sources. 



Linear STEP fit 




-0.10 - 
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Fig. 14. — Calibration of the RRG shear mea- 
surement method from simulated COSMOS im- 
ages containing a known input shear. The mea- 
sured shear on the y-axis includes the shear cali- 
bration factor C. Squares show measurements of 
7i , diamonds show measurements of 72 . The solid 
line is a linear fit to deviations from the ideal case 
of Tmeasured = 7in P ut for all points. The dashed 
line is a quadratic fit demonstrating that the cur- 
vature terms are negligible. 

The recovered shear measurement from the 
simulated data is presented in Figure Q3] We 
find that, in order to correctly measure the in- 
put shear on COSMOS-like images, the RRG 
method requires an overall calibration factor of 
C= (0.86™)-\ so that 



(14) 



The necessity for such calibratio n factors has 
long been known in th e field (e.g. iBacon et al 



2001uErben et al.ll2001r ). and is in accord with re - 
sults from STEP1. STEP2 and lHighet all (|2007l ) 
suggest that this may be intrinsic with KSB- 
related methods and, furthermore, that the cali- 



(c) 


( 2.1 ± 40.0) x 10" 4 


Cl 


( 5.6 ± 28.5) x 10- 4 


C2 


(-1.3 ± 28.0) x 10~ 4 


(m) 


(-0.3 ± 9.4) x 10~ 2 


m 1 


(-12.2 ± 6.5) x 10~ 2 


m 2 


( 11.6 ± 6.3) x 10" 2 



Quadratic STEP fit 

(c) (-23.3 ± 121.8) x 10~ 4 
(m) ( 20.1 ± 51.5) x 10~ 2 
jq) -2.84 ± 4.59 

Table 5: — Calibration of the RRG shear mea- 
surement method on simulated COSMOS images 
containing a known shear, described using STEP 
parameters. Figures are supplied after the appli- 
cation of the shear calibration factor. 



bration can vary for the two components of shear. 
For this reason, we fit each component separately, 
and use final calibration factors of C\ = (0.80) -1 
for 71 and C 2 = (0.92) _1 for 72. After this re- 
calibration, Table [5] shows STEP-likc estimates of 
the additive bias (c) and the multiplicative bias 
(m) obtained by fitting deviations of the recov- 
ered shear from the input shear. Both of these 
are consistent with ideal shear recovery, although 
the error on these estimates will be propagated 
through subsequent analyses. 

For the kind of 3D shear analysis for which 
COSMOS is so well-suited, the simultaneous cal- 
ibration of shears from an entire population of 
gala xies is insuffic i ent. S ince our companion pa- 
per, 



Massev et al 



(|2007j ). is concerned with the 
growth of the shear signal as a function of red- 
shift, it is crucial that the shear calibration be 
equally precise for both distant and relatively 
nearby galaxies. Given that more distant galax- 
ies are fainter and smaller and that the details of 
the shear measurement depend upon a fixed PSF 
size, pixel size, and noise level, this requirement 
is not trivial. We have therefore split the sim- 
ulated galaxy catalog in half by magnitude (at 
F8UW = 25.4) and by size (at d = 5.0 pix- 
els), and repeated the analysis. We find that our 
shear calibration (m) is robust for galaxies of dif- 
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ferent fluxes within 1% and of different sizes within 
4%. Redshifts were not available for the simulated 
galaxies, so a direct split in redshift was not pos- 
sible. Although this effect is clearly small in the 
regime of our current measurements, it will be sig- 
nificant in future weak lensing surveys where the 
error budget will be dominated by systematic un- 
certainties. 

5.5. Error on the Shear Estimator 

For each galaxy, the error on the measured 
shear is estimated using the same method as im- 
plem ented in the Photo pipeline i Lupton et al.l 
120011 ) to analyze data from the Sloan Digital Sky 
Survey. For each object, we assume that the opti- 
mal moments are the same as the moments corre- 
sponding to a best fit Gaussian. This formulation 
allows us to determine the covariance matrix of 
the moments in the usual way of non-linear least 
squares. To be precise, we model the image with 
a 2D elliptical Gaussian model 



M(x) = 



2tt|Q|V2 



exp[--(x-/i) T Q -/")]■ 



(15) 



This model has six parameters: the flux, /, the 
two centroids, /z, and the three moments which 
form the elements of the symmetric matrix Q. 
These parameters are noted pi. Next, we derive 
the \ 2 m t ne usual way 



X 2 ^EW x ij)-^] 2 , (16) 



where a is the sky noise level. We can now com- 
pute the 6x6 Fisher matrix which is the matrix 
of the second derivatives of the \ 2 



Fi 



ki 



1 d 2 x 2 



2 dpkdpi 

1 ^ <9M( Xij ) <9M(x y ) 



cr 2 4^ dp k 



dpi 



dp k dpi 

As is customary, we drop the second term which is 
proportional to the residuals. This term is usually 
very small compared to the first. The covariance 



matrix of the Gaussian parameters are then the 
inverse of the Fisher matrix. The 3x3 block of 
the covariance matrix corresponding to the sec- 
ond moments can then be extracted. Because the 
whole 6x6 matrix was inverted, this correctly 
marginalizes over centroid errors and other model 
parameter degeneracies. In this way, it differs from 
formulas which assume a constant (non-adaptive) 
weighting function and perfect centroiding. 

This Fisher matrix does not depend explic- 
itly on the data; it only depends on the best 
fit parameters and can be computed analytically. 
It is simply a function of four numbers: the 
three second order moments defined in Equation [7] 
and the signal-to- noise ratio (/ /er), which can be 
parametrized by the magnitude error (for exam- 
ple, magerr_AUTO) . Since the flux computed 
with SExtractor is not exactly equal to / (which 
would be the best fit Gaussian amplitude) we al- 
low for a single calibration factor and multiply the 
covariance matrix by this. We calibrate this factor 
with image simulations and verify that the errors 
are correctly predicted. 

Since the ellipticity components are computed 
from the moments, the variances of the ellipticity 
components can be computed by linearly propa- 
gating the covariance matrix of the moments. Fi- 
nally, the two ellipticity components can be shown 
to be uncorrelated with each other. 

6. Final Galaxy Selection 
6.1. Lensing Cuts 

Estimations of the gravitational shear will be 
improved by averaging only those galaxies with 
precise shape measurements on the condition that 
no ellipticity selection bias is introduced by the 
"lensing cuts" . We apply strict cuts to the Ci cat- 
alog that are designed to extract a sample of re- 
solved galaxies with reliable shape measurements. 
The resulting catalog is referred to as C4. Our 
"lensing cuts" are summarized in Table [6] and are 
based on the four following parameters: 

1. The estimated significance of each galaxy de- 
tection, where the significance is defined as 
s/n = flux_auto/fluxerr_auto, 



2. The first order moments, I xx and I., 



yy 



3. The total ellipticity, e = yef + e 
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4. The galaxy size as defined by the RRG d 
parameter (see §5.2[) . 

The final size cut is designed to select galaxies 
with well resolved shapes. Indeed, PSF correc- 
tions become increasingly significant as the size of 
a galaxy approaches that of the PSF and the in- 
trinsic shape of a galaxy becomes more difficult to 
measure. In COSMOS images, the typical size (as 
defined by d) of a star is about d* — 2.2 pixels 
(0.066"). Our size cut is thus equivalent to select- 
ing galaxies with d 9 > 1.6 x d-k. Note that in this 
section, as well as in all following sections, d has 
not been corrected for the PSF. 

The ellipticity cut at e < 2 may be surprising 
given that, by definition, cllipticities are restricted 
to e < 1. In reality, however, because of noise, 
it is possible to measure an ellipticity of e > 1. 
Selecting galaxies with e < 1 could introduce an 
unwanted ellipticity bias, but, because we are only 
interested in ensemble averages, an acceptable so- 
lution is to cut out only a small number of large 
outliers (e > 2). 

Note that further cuts may be required for 
some applications (for example, to isolate only 
those objects with well-measured photometric red- 
shifts). In addition to these cuts, a galaxy- by- 
galaxy weighting scheme may also be used to min- 
imize the impact of shape measurement noise. 

6.2. Effective Galaxy Number Density 

Once the "lensing cuts" have been applied, 
the C4 catalog contains only those galaxies which 
are useful for lensing analyses with the COSMOS 
data. Predictions for future weak lensing surveys 
based on COSMOS results must consider the "ef- 
fective number" of galaxies that are actually useful 
for lensing purposes (and not the number of raw 
detections). With these considerations in mind, 
we define the "effective galaxy number density" , 
N g (z), as the total number of galaxies within C4, 
per unit area, and with redshifts below a given red- 
shift, z. The equivalent quantities as a function of 
galaxy magnitude and size are N g (m) and N g (d). 
For each galaxy, the SExtractor MAG_AUTO pa- 
rameter is used to estimate the magnitude and the 
RRG d parameter is used as an estimate of the 
size. We also consider the derivatives of N g (z), 
N g (m), and N g (d) so that: 




N g {m) = / n g (m')dm' , (18) 
J20 



N g (d) = / n g (d')dd'. (19) 

Jd 

The total number density of galaxies in C4 is 
noted as N g and is significantly lower than the 
number density of detected galaxies. In total, the 
final lensing catalog C4 contains 3.9xl0 5 galaxies 
with accurate shape measurements and 2.8xl0 5 
galaxies with both shape and photometric red- 
shift measurements. Table [5] shows a summary 
of the different steps leading to this catalog. The 
surveyed area of COSMOS is 1.64 deg 2 leading 
to an overall number density of N g ^66 galaxies 
per arcmin 2 for the first sample and ~50 for the 
second. These numbers can be contrasted to a 
more sparse 15-25 galaxies per arcmin 2 typically 
resolved with deep, ground-based surveys. In Fig- 
ures [15] and [16] we show the effective densities de- 
fined above as well as their corresponding deriva- 
tives. From these figures, we can draw the follow- 
ing conclusions: 

• Over 60% of the COSMOS source galaxies 
are at redshifts higher than z = 1. The 
COSMOS weak lensing data is therefore a 
powerful probe of the dark matter distribu- 
tion from z ~ 1 to the present day. 

• About 18 galaxies per arcmin -2 (73%) are 
discarded when we select only those galax- 
ies with accurate photometric redshifts. The 
primary cause of this loss are the larger ar- 
eas masked out in the ground-based data as 
compared to the ACS data. Future space- 
based weak lensing surveys could recover the 
remaining 27% by using space-based, multi- 
wavelength imaging to derive photometric 
redshifts. 

• The effective number density rises very 
steeply with decreasing galaxy size; over 
50% of our total number of sources have 
d < 5 pixels (0.15"). The current size cut 
for the COSMOS data is d = 3.6 pixels and 
the size of the ACS PSF is d* = 2.2 pix- 
els (0.066"). By pushing this size barrier 
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Table 6 
Lensing cuts applied to C 2 



Parameter 


Galaxies retained in C4 


I xx and I yy 


Finite* 1 I xx and I yy 


RRG size parameter d 


d > 3.6 pixels 


Significance 


S/N > 4.5 


Total cllipticity c e 


e < 2 



indicating that the RRG code converged 
b Uncorrected for the PSF 
"Corrected for the PSF 



to even smaller values, future surveys could 
very quickly obtain much higher effective 
densities. 

Finally, the cuts that we have applied to the 
COSMOS data (in particular the size cut) are 
stricter than would be necessary without signifi- 
cant CTE effects that degrade the shape measure- 
ments of the faintest galaxies. Next generation 
space-based missions designed to avoid the prob- 
lems encountered and identified in COSMOS, as 
well as future implementations of the COSMOS 
catalog, will undoubtedly achieve higher effective 
number densities. 

7. Intrinsic Shape Noise: a Fundamental 
Limit to Weak Lensing Measurements 

Under the assumption of weak gravitational 
lensing, a source galaxy with intrinsic shape Si n t 
and observed ellipticity e a bs is related to the grav- 
itational lensing induced shear 7 according to: 

Sobs = Sint + 7- (20) 

Throughout this paper, the gravitational shear 
is noted as 7 whereas 7 represents our estimator 
of 7. The above relationship indicates that galax- 
ies would be ideal tracers of the distortions caused 
by gravitational lensing if the intrinsic shape £{ n t 
of each source galaxy was known a priory. A 
quick glance at an ACS image however, reveals 
that galaxies display a very wide variety of shapes 
which unfortunately prevents the extraction of 
7 for any single galaxy. Lensing measurements 
thus exhibit an intrinsic limitation, encoded in the 
width of the ellipticity distribution of the galaxy 
population, noted here as o~i n t, and often referred 



to as the "intrinsic shape noise". Because the 
shape noise (of order <jj nt ~ 0.26) is significantly 
larger than weak shear (typically 7 ~ 0.01 for cos- 
mic shear) , 7 must be estimated by averaging over 
a large number of galaxies. In this case equation 
f2ni simplifies to: 

< e obs >=< 7 > . (21) 

The uncertainty in the shear estimator, 07, 
arises from a combination of unavoidable intrinsic 
shape noise, af nt =< sf nt > and the measurement 
error of galaxy shapes 0-^ eas : 

a \ = a int + a meas- ( 22 ) 

In the following analysis, 07 will be referred to 
as the shape noise and o~i n t will be called the in- 
trinsic shape noise. The former includes the shape 
measurement error, <r meas , and hence will vary ac- 
cording to the data-set as well as the shape mea- 
surement method that is employed. Note that the 
uncertainty contributions from, photon noise, PSF 
correction, CTE calibration, and the shape mea- 
surement method are all included in our definition 
of cr„ ea;s . The weak lensing distortions averaged 
over the whole COSMOS field arc small, and rep- 
resent a negligible perturbation to equation 1221 

Instead of the simple arithmetic mean of equa- 
tion [21] many lensing practitioners adopt in some 
form or another, an optimized weighting scheme in 
order to estimate 7 which often incorporates both 
the measurement error and t he shape noise (for 
e.g., iBernstein k. Jarvisl l2002h . A constant value 
(of order 0.3) is often assumed for cr int . However, 
it would not be surprising that the same processes 
that shape galaxy formation also lead to a vari- 
ation of a int as a function of magnitude, galaxy 
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Fig. 15. — In the upper panel we show the ef- 
fective number density of galaxies as a function of 
F81AW magnitude. Only resolved galaxies with 
precise shape measurement are included in these 
counts. In the lower panel, we show the effec- 
tive number density of galaxies as a function of 
redshift. The effective density still evolves sharply 
after z > 1 demonstrating that the COSMOS lens- 
ing data is a powerful probe of structures at z < 1. 

type, or redshift. Furthermore, as weak lensing 
surveys increase in both scale and depth, it of in- 
tense interest to obtain accurate estimates of the 
intrinsic shape noise floors that these surveys must 
confront. For these reasons, we undertake a mea- 
surement of the shape noise, a^, as well as the 
intrinsic shape noise, <7j n t, as a function of magni- 
tude, size and redshift. A more detailed analysis 
of the intrinsic shape noise as a function of galaxy 
morphology will be the subject of a future paper. 

First, we estimate the shape noise <x^ directly 
from the COSMOS data as a function of size and 
F814VF magnitude. To derive <jj we consider: 



Fig. 16. — In the upper panel we show the ef- 
fective number density of galaxies as a function 
of size. The dotted vertical line indicates the 
size cut that we make in order to extract galax- 
ies with precise shape measurements. This size 
cut is such that d g > 1.6 x d-k. N g rises very 
sharply as a function of decreasing d demonstrat- 
ing that small galaxies make up the majority of 
our lensing sources. The lower panel shows N g for 
COSMOS galaxies with accurate photometric red- 
shifts and for magnitude cuts of F814VF < 25 and 
F8UW < 24.5. 

• The mean variance of both shear compo- 
nents (including correction factors), a-y = 
(ayi + cr^ 2 )/2 

• Galaxies with well measured shapes (see § [6]) 

Second, we derive an empirical estimation of the 
shape measurement error, a meas , using a sample 
of 27000 galaxies that belong to overlapping re- 
gions of adjacent pointings. Each of these galax- 
ies provide us with two independent shape mea- 
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Fig. 17. — The observed scatter in the shear as a 
function of F814W magnitude and size. The scat- 
ter is a combination of intrinsic shape noise, cr inty 
and a shape measurement error, a meas - The shape 
measurement error is determined by a theoretical 
model and tested using galaxies from overlapping 
regions. Shaded regions indicate the RMS width of 
the measurement error distribution. The intrinsic 
shape noise appears to increase slightly as a func- 
tion of magnitude but is independent of size. At 
faint magnitudes and small galaxy sizes, the shear 
scatter increases rapidly due to large measurement 
errors. The dashed vertical line indicates the size 
and magnitude cut that we apply in order to con- 
struct Figure [T51 



surements. Using these overlaps, we find that the 
shape measurement error is a function of both size 
and magnitude and increases beyond a mea s = 0.1 
for MAG_AUTO > 24.5 and d < 7. Third, we com- 
pare this empirical determination of cr meas to the 
theoretical one derived in § 15.51 We find that the 
theoretical model of § 15.51 does remarkably well in 
predicting <7 meas as a function of size and magni- 
tude. Thus confident in the validity of this model, 
we adopt it for subsequent derivations. Finally, 
using equation [22] and the shear measurement er- 
ror a meas , we extract the intrinsic shape noise of 
our galaxy sample as a function of size, magnitude 
and redshift. The results are shown in Figures IT71 
and [13 

As can be seen in Figure I17i large measure- 
ment errors lead to an increase of at small sizes 
and faint magnitudes. The intrinsic shape noise 
however, appears to change little with either size 
of magnitude and remains constant at a value of 
&int ~ 0.26. The slight apparent increase of the 
intrinsic shape noise at fainter magnitudes is prob- 
ably due to the simplified measurement error esti- 
mator that we are using (indeed, the overlaps in- 
dicate slightly higher errors). From this analysis, 
we can draw the following conclusions: 

• The intrinsic shape noise varies little from 
z = to z = 3. Deep space based weak lens- 
ing surveys will therefore confront equivalent 
intrinsic shape noise floors as their shallower 
counterparts. 

• Measurement errors lead to an increase in 
the shape noise as a function of size and 
magnitude. A joint improvement in both 
imaging quality as well as shape measure- 
ment methodology will lead to shape noises 
that are closer to the intrinsic floor of 0.26. 

• We have yet to explore if the intrinsic shape 
noise varies as a function of galaxy type. If 
so, shear measurements could be improved 
by incorporating a galaxy-type discriminant 
into the weighting scheme. This will be the 
focus of a future paper. 

8. Conclusion 

We have carefully constructed a weak lens- 
ing catalog from 575 ACS/HST tiles, the largest 
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Fig. 18. — Intrinsic shape noise as a function 
of photometric redshift. Galaxies have first been 
selected to have measurement errors less than 0.1 
(d > 7 and MAG_AUTO < 24.5) so that the scatter 
in the observed shear closely matches the intrinsic 
shape noise. A linear fit to the shape noise as a 
function of redshift reveals a flat distribution with 
a mean of ai n t — 0.26. 

space-based survey to date. We have established 
the quality of this catalog and analyzed the COS- 
MOS redshift distribution, showing broad agree- 
ment with other photometric redshift surveys out 
to z ~ 1. The photometric redshifts are currently 
limited by the lack of deep if-band data but 
will rapidly improve as this data soon becomes 
available. Shapes have been measured for over 
3.9 x 10 5 galaxies and corrected for distortions in- 
duced by PSF and CTE effects. Simulations have 
been used in order to calibrate our shear measure- 
ment method and a STEP-like analysis has been 
performed, demonstrating our ability accurately 
measure shear with negligible additive and multi- 
plicative bias. The effective number density of the 
COSMOS weak lensing catalog is 66 galaxies per 
arcmin~ 2 (48 when we consider only those with 
accurate photometric redshifts). A large fraction 
of these galaxies are at z > 1 making COSMOS a 
powerful probe of the dark matter distribution and 
its evolution from z = 1 to the present day. The 
COSMOS survey is also of foremost importance 
for the preparation and design of future wide field 
space-based lensing missions. Regarding the de- 
sign of such missions, our main conclusions from 
working with the COSMOS data are the following: 



1. Understanding and correcting for the time 
varying PSF and calibrating CTE effects 
were two of the most difficult challenges en- 
countered with the COSMOS data. Reduc- 
ing these two systematic effects should be 
a key specification in the design of next- 
generation telescopes and instruments. 

2. Because 1) small galaxies are not as read- 
ily detected from the ground and 2) because 
larger areas are masked out in the ground- 
based data than in the ACS imaging, we lose 
27% of our source sample when we select 
only those with accurate photometric red- 
shifts. 

3. The effective number density of galaxies is a 
very sensitive function of survey depth and 
resolution. The capability to resolve and ac- 
curately measure the shapes of very small, 
faint galaxies will be key in obtaining num- 
ber densities of over 66 galaxies arcmin . 

4. Finally, we have derived the intrinsic shape 
noise of the galaxy sample and demonstrated 
that it remains fairly constant as a func- 
tion of size, magnitude and redshift. Weak 
lensing measurements with deep space-based 
imaging are therefore on par with more shal- 
low imaging in terms of the intrinsic shape 
noise floors that they must overcome. 

The COSMOS weak lensing data described 
in this paper has already been used to mea- 
sure cosmological parameters and to demonstrate 
the feasiblity of t he "tomography" technique 
(jMassev et al.l 120071 ). Striking weak lensing mass 
maps of the COSMOS field have been made that 
reveal tantalizing evidence of a complex interplay 
between the baryon and the dark matter distri- 



bution (jMassev et alJ l2007ah . Yet more lensing 



analyses are underway including a galaxy- galaxy 
lensing and a group-galaxy lensing study that will 
undoubtedly lead to a better understanding of the 
relationship between baryonic and dark matter 
structures and of its evolution over cosmic time. 
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